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Abstract 

In this note we aim to infer a model for the response of a Pierre Auger water-Cherenkov 
detector to an ideal single muon. The main goal of this analysis is to provide analytical support 
for muon counting techniques. In this note we derive the probability distribution of the muon 
tracklength as a function of the zenith angle of the muon. 



1 Introduction 


This note is the first in a series of nofes in which we aim fo infer a model for fhe response of a 
Pierre Auger wafer Cherenkov defecfor {tank from now on) fo an ideal single muon. The single 
muon response is a subjecf fhaf was fhoroughly explored in fhe early phase of fhe Auger col- 
laborafion [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. The main purpose of fhese sfudies 
was fo undersfand fhe mean muon response (more precisely, fhe so-called muonic peak or fhe 
"muon hump") in order fo define an SD energy esfimafe. Firsf, based on fhe muonic peak model, 
a calibration procedure was designed for estimating fhe fofal signal in individual fanks. Then fhe 
fofal signals were combined fo compufe one observable per shower, which was finally calibrafed 
fo fhe energy esfimafe of fhe fluorescence defecfor (FD). 

Since fhe SD energy esfimafe is calibrafed fo fhe FD energy esfimafe, fhe procedure is rela¬ 
tively robusf againsf sysfemafic biases of fhe mean muon response esfimafe as long as fhe biases 
are sfable among differenf fanks and nof changing wifh fime. The flucfuafions of fhe muonic 
signal are also nof foo imporfanf from fhe poinf of view of fhe energy esfimafion. On fhe ofher 
hand, fhe goal of muon-counfing [17, 18,19, 20] is fo design a muon densify esfimafor wifhouf fhe 
need for oufside calibrafion. Now biases due fo fhe fank geomefry and fhe energy-dependence 
of fhe muonic response become imporfanf. Furfhermore, since fhese fechniques are based on 
cerfain sfafisfics of fhe individual muonic "jumps", undersfanding fhe flucfuafions of fhe muon 
response is of great importance. Our ultimate goal is similar to the program outlined in [12]: ob¬ 
tain a full paramefrizafion of fhe muonic signal fhaf can be used in a Monfe Carlo Markov chain 
[21] reconsfrucfion approach as well as for fine-funing fhe muon counfing fechniques. Beside 
fhis principal goal, we believe fhaf fhis refined model may also help fo improve fhe SD energy 
esfimafe (mainly by decreasing fhe sfafisfical error on fhe individual shower esfimafes along fhe 
lines of [22]). The obfained model may also serve as a basis for a toy Monte-Carlo tank simulator 
that can be used to quickly generate a large number of fank signals. 

1.1 Introduction to this note 

In a crude model, the total muonic signal is proportional to the tracklength of the muon in the 
water tank, so the fluctuation of the tracklength at a given zenith angle is the main source of 
the fluctuation of the total muonic signal. In this note we derive the tracklength distribution 
I 6) of a muon that crosses the tank at zenith angle 6. Most of the formulas can also be 
found in [11]. Our main contribution is that we give clean formulas that can be used directly in 
a probabilistic generative model. We also add one term for muons that enter and leave at the 
lateral of the tank. This part of the signal, which was omitted by [11], may be important for 
inclined showers. 


2 The distribution of the muon tracklength 


In this section we derive the tracklength distribution 'Pi{L \ 0) of a muon that crosses the tank 
zenith angle 0. To transform the quantities into unit-less variables, we first apply the substitutions 


i = 


L 

2R 


and 



( 1 ) 


where h = 1/3 for an Auger tank with height H = 1.2 m and radius R = 1.8 m. From now on 
we will only derive the distribution 7?(^ | 0) as a function of the unit-less variable £. The original 
distribution in L can be recovered by the transformation 

Vdl I 0) = 4p(L/2R I 0). (2) 


The distribution can be decomposed into three components: the lid-base (lb) term of the form 
Vlhd I ^(lt> I 0) representing muons entering at the lid of the tank and leaving at the base, the 
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lid-lateral (li) term | 9) V{\i \ 6) representing muons entering at the lid and leaving at the 
lateral (or, symmetrically, entering at the lateral and leaving at the base), and the lateral-lateral (la) 
term \ 9) V{\a \ 9) representing muons entering and leaving at the lateral (for these regions 
see projected tank schematic in Fig. 1). We have thus 

I 9) = I 9) V{lh I 9) + riiH \ 9)V{li I 9) + Pi,(^ I 9) P(la I 9). (3) 



Figure 1: Projection of fhe Auger fank along fhe incoming muon direcfion af F = 50°. The fhree 
componenfs of fhe fracklengfh disfribufion p(^ | 50°) are visible: fhe lid-base (lb), fhe lid-lateral 
(li), and fhe lateral-lateral (la) ferms. 


2.1 Case probabilities 

We proceed by firsf compufing fhe case probabilities VlVo \ 9), T’(li | 9), and T’(la | 9). First we 
need the area A{9) of fhe fank projecfion for zenifh angle 9, 

A{9) = nR^\cos9\+2RHsin9. (4) 


Nofe fhaf fhe absolufe value is used in order fo render fhe equafions useful also for cases when 
muons are enfering fhe fank wifh zenifh angle 9 larger fhan 90° (i.e., useful for albedo sfudies). 
Measuring area in unifs of fhe fop surface At = ttR^ and using fhe unit-less h, we can define 


a{9) 


m 


cos0| -H ^/zsinF 


||cosF|(f -F2 m), 


(5) 


where 

u = h\ fanF|. (6) 

Measuring fhe partial areas (see Fig. 3(a)) in the same relative unit Axx = ^xx/ At, we obtain for 
fhe lid-base parf 


flib(F) 


^1 cosF| arccos(/zfanF) — hsm9y/l — (/zfanF)^^ 

0 


i arccos u — uy/l — u?- 
0 


if M ^ 1, 
ofherwise. 


if I fanF| ^1/h, 
ofherwise. 


(7) 

( 8 ) 
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for the lid-lateral part 


%(0) =2(1 cos0| -flib(0)) = 


= ||cos0| 




arcsin u 


+ w\/l — if w ^ 1, 


TT 

2 


ofherwise. 


and for fhe laferal-laferal parf 

“la(^) = ^hsinO- |cos0| +flib(0) = 


2u — arcsin w — u^/l — u?- if m ^ 1, 


= ||COS 0 | 


2w - 


ofherwise. 


The case probabilities (see Fig. 3(b)) are fhen given by 


V{\h I e) = 


p(ii|0) = ^, and ^(la I 0) = 


a{0) ' \ \ 

Nofe fhat for fhe final case probabilities, fhe 2| cos 0| /tt ferm cancels in fhis process. 


(9) 

( 10 ) 


( 11 ) 

( 12 ) 


(13) 


2.2 Probability distribution functions 

The lid-base ferm Pib(^ | 0) is a simple Dirac delfa 


V\h{^ I - ^max). 


where fmax(0) is the maximum tracklength (see Fig. 2) defined as 




fi/|cos0| if I tan0| ^ 1//!, 
1/ sin0 otherwise. 


(14) 


(15) 


Note that the case 1 / sin 0 is never in effect in the lid-base term since the case probability 
P(lb|0) is zero in this range. The lid-lateral and lateral-lateral terms are 


Pli(^ I 0) 


Pla(^ I S) 


1 

flli(0) 


f sin20| v^l — (^sin0)2 

[o 


ifO^^^ W(0), 

otherwise. 


1 

«la(0) 


4 .fsin^ 6{h — ^\ cos0|) 
" y/l — (£ sin 0)2 

0 


ifO^£^^inax(0), 

otherwise. 


(16) 


(17) 


respectively. 
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Figure 2: The mean tracklength L{6) = E {'pi{L\9)} (red), the maximum tracklength Linax(^) = 
2K^max(F) (green), and the region that contains 70% of the probability mass (blue) as a function 
of fhe zenifh angle 6. 


2.3 Cumulative distribution functions 

The cumulative distribution functions of the continuous terms are 
■ I cos 0\ 


[ I o)de' 

Jo 


^ aii(F) 


arcsin(t' sin 9) + £ sinF-y/l — (^sinF)^ 
arcsin v + v\/l — 


arcsin u + mVI — 
arcsin u + uVl — 


71/2 


if u ^ 1, 
otherwise. 


/Vla(^' I 0)d£' 

Jo 


I cos 9\ 

flla(0) 


2ti| tan0| — arcsin(^sin0) — {2h\ tanF| — tsinF) \J'^~ (tsinF)^ 
2w — arcsin u — (2w — v)Vl 


where 


2w — arcsin M — wVl — 

2u — arcsin u — (2w — v)^/\ —v^ 
2u — 7t/2 

u = ^ sin 9, 


if u ^ 1, 


otherwise. 


(18) 

(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 
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and the full cumulative distribution is 



arcsin v 


+ 2m + (3u — 2m)\/ 1 — i 


2m + 7r/2 


if u ^ min(l,M), 
otherwise. 


( 24 ) 




e n e [“] 

(a) relative areas (b) case probabilities 

Figure 3: (a) The relative partial areas Mib(F)/ ^iii(ll)/ ^iia(^)/ ^^id their sum a{6), and (b) the case 
probabilities ViVo \ 9), V{li \ 9), and V{la | 0) as a function of the zenith angle 9. 


2.4 Mean tracklength 

It has already been noted several times that the mean tracklength L(F) (see Fig. 2) for muons 
arriving with a specific zenith angle 9 is quite easy to obtain without integrating L'pi{L \ 9). If 
we enlarge the puncturing track of a muon so that it has the same cross section dA along the 
path, a small volume dV is obtained. Adding up all such parallel tracks will eventually amount 
to the whole volume of the tank, 

V = JdV = JldA. (25) 


The volume of the whole tank is V = nRp-H and the right integral of Eq. (25) is in fact equal to 
the average tracklength multiplied by the projected area f dA = A(9). Therefore V = L(9)A(9) 
and 


L(9) = E 9)} = 


V 


H 


2Rh 


I cos 01 + ^ sin 9 


I cos 9\ + ^ sinF 


(26) 


3 Conclusion 

In this note we derived the probability distribution of the muon tracklength as a function of its 
zenith angle 9. The results show that the mean tracklength is relatively stable in the [0°,60°] 
interval, however the fluctuations increase rapidly for 9 > 10° (see Fig. 4). For inclined showers 
{9 > 60°), the mean tracklength also increases by a factor of 2. 
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Figure 4: The relative standard deviation ^/Var {'pi(L \ 9)}/L{6) of the tracklength as a function 
of fhe zenith angle 6. 
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